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ABSTRACT 



Context. In this paper we analyse the spectra of D-type SS Hl-36 within a colliding-wind scenario. 

Aims. We aim to analyse the properties of this object taking into account the observational data along the whole electromagnetic 
spectrum, in order to derive a self-consistent picture able to interpret the nature of the system as a whole. 

Methods. After constraining the relative physical conditions by modelling more than 40 emission lines from radio to UV, we are able 
to explain the continuum spectral energy distribution by taking into account all the emitting contributions arising from both the stars, 
the dust shells and the gaseous nebula?. 

Results. A comprehensive model of the radio spectra allows to reproduce the different slopes of the radio profile and the turnover 
frequency, as well as the different size of the observed shocked envelope at different frequencies in the light of the different contri- 
butions from the expanding and reverse nebulae. The IR continuum unveils the presence of two dust shells with different radii and 
temperatures, which might be a distinctive feature of D-type symbiotic systems as a class of objects. The broad profiles of IR lines 
direct us to investigate whether an X-ray jet may be present. 

Conclusions. This insight leads us to indicate Hl-36 as a promising X-ray target and to encourage observations and studies which 
consistently take into account the complex nature of symbiotic stars throughout the whole electromagnetic spectrum. 

Key words, binaries: symbiotic - stars: individual: Hl-36 



1. Introduction 

Symbiotic systems (SS) are currently understood as interacting 
binaries composed of a compact star, generally but not neces- 
sarily a white dwarf (WD) which is the source of the ionizing 
radiation; a cool giant star, which is at the origin of dust forma- 
tion and ejection episodes; and different emitting gas and dust 
nebulae. The large variety of the properties of these systems of- 
ten makes their classification as a single group difficult; however, 
these differences are probably more related to distinct physical 
conditions like masses, orbital periods and temperature of the 
hot source, than different kind of phenomena. 

Because of the associated circumstellar matter and energetic 
activity, many SS emit detectable radiation across nearly the 
entire electromagnetic spectrum. On the basis of near-infrared 
(NIR) colours, SS were classified in S and D types (Webster & 
Allen 1975) according to whether the cool star (S-type) or dust 
(D-type) dominates the 1-4 fim spectral range. As with regards 
to the continuum radio emission, it has been modelled as ther- 
mal bremsstrahlung radiation arising from (a) the photoionized 
component of the cool star wind (Seaquist et al. 1984, hereafter 
STB - Taylor & Seaquist 1984), (b) a wind associated with a 
slow nova-like eruption from the hot companion, (c) the interac- 
tion region involving both winds in mutual collisions (Watson et 
al. 2000, Kenny & Taylor 2005, Bisikalo et al. 2006). Whatever 
the case, it is widely recognised that the radio spectra provide 
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unique informations about the mass-loss phenomena at the basis 
of symbiotic activity. Eventually, there is growing evidence that 
SS can be the site for jet structures and synchrotron emission, 
confirming the primary role of shocks in the nowadays widely 
accepted interpretation of SS as colliding-wind binary systems. 

Hl-36 was among the most interesting objects in the exten- 
sive Purton's radio survey of emission-line stars (1982). The par- 
ticular attention that was drawn to it since the Haro's discovery 
(1952) was due to a curious set of misunderstood and genuine in- 
sights: catalogued as planetary nebula despite its extremely high- 
excitation emission line spectrum and its imposing infrared ex- 
cess, Hl-36 was even thought to be, now we know erroneously, 
the optical counterpart of the Uhuru hard X-ray source 3U 1746- 
37 (Giacconi et al. 1974). Finally, Allen permanently classified 
Hl-36 as a D-type SS by including it in his catalogue (Allen 
1981) and dedicating to it a paper (Allen 1983) that so far still 
represents the most complete observational work on such object. 

One of the most important results emerging from that study 
is that the Mira star experiences an extraordinary reddening 
(A,, ~ 20 mag.), especially when compared to the one towards 
the emission-line regions (only 2.2 mag.): this allowed Allen to 
point out that the cool component is heavily embedded in its 
own circumstellar dust shell, illuminated from the outside by a 
T* ~ 150.000 K star. The derived cool star and dust temperatures 
would then be, respectively, T M ~ 2500 K and T d ~ 700-800 K. 
Another interesting feature of Hl-36 lies in the radio range: as a 
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matter of fact, it is one of a small number of radio sources whose 
spectra are fiat at high frequencies but turn towards a spectral 
index near +1 at lower frequencies (Purton 1977). Several 
papers hence tried to explain this behaviour in terms of free-free 
emission from the ionized cool star wind, unfortunately without 
taking into account the possibility of any outflow from the 
compact star (for a further description of this models, see Sect. 
4.2.1). 

The following papers, which presented both imaging (e.g. Bhatt 
& Sagar 1991; Corradi et al. 1999, who resolved the complex 
nebula at radio and optical wavelengths) and new spectra (e.g. 
Costa & de Freitas Pacheco 1994; Pereira 1995; Pereira et 
al. 1998), have not substantially been modifying the original 
scenario drawn by Allen's analysis, whose physical parameters 
have still to be considered as the most reliable ones available in 
literature. Interestingly, Hl-36 is the only SS known to support 
OH- (Ivison et al. 1995), as well as SiO- (Allen 1989) and H 2 0- 
masers (Ivison et al. 1998). 

In this paper we aim to model Hl-36 (Tabled by combining 
literature observations from radio to UV collected over a period 
of about 30 years (Sect. 2). Within a colliding-wind scenario and 
by means of the SUMA code (Sect. 3), we start modelling more 
than 40 emission lines from IR to UV in order to constrain the 
physical conditions across the whole system (shocked nebulae, 
dust shells, hot and Mira stars). In the light of the derived param- 
eters, we are then able to fit the composite continuum SED in a 
self-consistent way. The results of this cross-checked method are 
extensively presented in Sect. 4. Concluding remarks appear in 
Sect. 5. 

2. Observational details 

2.1. Radio-mm observations 

In Table|2]we present the data we collected from literature in or- 
der to ensure as much a complete spectral coverage in the radio- 
mm range as possible. 

The most important references about the Hl-36 radio spectrum 
remain Purton's papers (Purton 1977, 1982), in which observa- 
tions from about 2 GHz to 90 GHz obtained with the CSIRO 
64-m telescope at Parkes, Australia, the NRAO 11-m telescope 
at Kitt Peak, USA, and the NRAO three-element interferome- 
ter at Green Bank, USA, are presented. In the next years several 
works, focused mainly on the mass loss properties of the cool 
stellar component (e.g. Jones 1985) and on the presumed corre- 
lations among different spectral bands (Seaquist et al. 1993), per- 
formed observations at intermediate frequencies. They allowed 
us to extract a well defined radio profile and eventually to con- 
firm that the object has not been significantly variable in the last 
30 years. 

This ensures that collecting data coming from a wide temporal 
range is, at least for this SS, a reliable approach. Summarising, 
we are able to present a radio-mm profile composed by 16 points 
from 843MHz to 230GHz (Fig. Q] bottom panel). 

2.2. Infrared data 

D-type SS are known to be variable systems in the IR range, 
where we observe the contributions of different emitting compo- 
nents, i.e. cool star and dust shells. Quite surprisingly, as in the 
radio case, the Hl-36 IR spectra do not show dramatic variations 
in the data flux over the years we are going to discuss (Table [3] 
Fig- ID top panel): indeed, IRAS points (1982) agree well with 



ISO-SWS spectrum (taken in 1996) and with the 34.6 yum mea- 
surement presented in He et al. (2005); also the accordance with 
the MSX6C photometric points (Egan et al. 2003) is reasonably 
good. At NIR wavelengths the Allen's NIR bands are slightly 
shifted with respect to the DENISE and 2MASS bands, espe- 
cially the H and J bands which feel the increasing contribution 
of the cool star. This is confidently due to the Mira variations, 
whose period is about 450-500 days (Whitelock 87). 

2.3. Optical spectra and UV spectrophotometry 

Also in this case the analysed data are a collection of different 
observing sessions from 1975 to 1979. These are presented in 
Allen (1983), where the reader is sent for further technical de- 
tails. Moreover, in that paper it is explicitly stated that the opti- 
cal data "are combined in the belief that variability of the lines 
is slight". This has been even more reconfirmed by recent spec- 
tra, which in fact did not show relevant changes in the line and 
continuum flux level (Costa & de Freitas Pacheco 1994, Pereira 
1995, Pereira et al. 1998). 
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Fig. 1. Observational data. Top panel: IR spectral range - for sake 
of clarity we don't show the ISO spectrum. Bottom panel: radio 
spectral range. 



3. The theoretical framework 

3.1. The colliding-wind scenario 

In the past years, theoretical models (Girard & Willson 1987, 
Kenny & Taylor 2005) as well as observations (Nussbaumer et 
al. 1995) have categorically shown that in SS both the hot and 
cool stars lose mass through strong stellar winds which collide 
within and outside the system, hence creating a complex net- 
work of wakes and shock fronts which result in a complicated 
structure of gas and dust nebulae (Nussbaumer 2000). 
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Ref. 



a (J2000) 
6 (J2000) 
WD temperature 
Cool star spectrum 
Mira temperature 
Mira pulsation period 
Dust temperature 

Distance 
Binary separation 
Orbital period 
Other names 



17 49 48.1 
-37 01 27.9 
1.5 10 5 K 
Mira - M4-M5 
2500 K 
450-500 d 
700 K 
4.5 kpc 
3 10 16 cm 
>10D 3/2 yr 
PK 353-04 l=Hen2-289= 
=Haro 2-36=IRAS 17463-3700 



Belczynski et al. (2000) 
Belczyfiski et al. (2000) 

Allen (1983) 
Allen 83 - Medina-Tanco & Steiner (1995) 

Allen (1983) 
Whitelock(1987) 

Allen (1983) 

Allen (1983) 

Allen (1983) 

Allen (1983) 
Belczynski et al. (2000) 



In this paper, as previously done for other SS (see Sect. 3.2), we 
refer to two main shocks: the shock between the stars facing the 
WD, which is a head-on shock (hereafter the reverse shock), and 
the head-on-back shock, which propagates outwards the system 
(hereafter the expanding shock). Both the nebulae downstream 
of the shock fronts are ionized and heated by the radiation flux 
from the hot star and by shocks. The photoionizing radiation 
flux reaches the very shock front of the reverse shock, while it 
reaches the edge opposite to the shock front downstream of the 
expanding shock. 

3.2. The calculation code 

The results presented in this work are performed by SUM A 
(Viegas & Contini 1994; Contini 1997), a numerical code 
that simulates the physical conditions of an emitting gaseous 
cloud under the coupled effect of ionization from an external 
radiation source and shocks, and in which both line and con- 
tinuum emission from gas are calculated consistently with dust 
reprocessed radiation (grain heating and sputtering processes 
are also included). The derived models have been successfully 
applied to several SS, e.g. AG Peg (Contini 1997, 2003), HM 
Sge (Formiggini, Contini & Leibowitz 1995), RR Tel (Contini 
& Formiggini 1999), He2-104 (Contini & Formiggini 2001), 
R Aqr (Contini & Formiggini 2003), HD330036 (Angeloni et 
al. 2007b), as well as to nova stars (V1974 Cyg, Contini et al. 
1997 - T Pyx, Contini & Prialnik 1997) and supernova remnants 
(Kepler's SNR, Contini 2004). 

The calculations start with gas and dust entering the shock 
front in a steady state regime: the gas is adiabatically compressed 
and thermalized throughout the shock front. In the downstream 
region the compression is derived by solving the Rankine- 
Hugoniot equations (Cox 1972): the downstream region is au- 
tomatically divided in plane parallel slabs in order to calculate 
as smoothly as possible the physical conditions throughout the 
nebula. Radiation transfer and optical depths of both continuum 
and lines are calculated for a steady state: in particular, radiation 
transfer of the diffused radiation is taken into account following 
Williams (1967). The fractional abundance of the ions in differ- 
ent ionization stages is calculated in each slab by solving the 
ionization equilibrium equations for the elements H, He, C, N, 
O, Ne, Mg, Si, S, CI, Ar, and Fe. The electron temperature in 
each slab is obtained from the energy equation when collisional 
processes prevail and by thermal balancing when radiation pro- 
cesses dominate. 

Compression downstream strongly affects the gas cooling rate 



by free-free, free-bound, and line emission: consequently, the 
emitting gas will have different physical conditions depending 
on the shock velocity and on the pre-shock density. For example, 
self-absorption of free-free radiation, consistently calculated in 
every slab, affects the emission spectral index at low radio fre- 
quencies (see Sect. 4.2.1) 

Dust is included in the calculations, too. Dust and gas are 
coupled throughout the shock-front and downstream by the mag- 
netic field. In each slab the sputtering of the grains is calculated, 
leading to grain sizes which depend on the shock velocity and on 
the gas density. The temperature of the grains, which depends on 
the grain radius, is calculated by radiation heating from the ex- 
ternal (primary) source and by diffuse (secondary) radiation, as 
well as by gas collisional heating. The dust reprocessed radia- 
tion flux is calculated by the Plank-averaged absorption coeffi- 
cient of dust in each slab, and integrated throughout the nebula 
downstream. 

The input parameters which characterise the shock are the 
shock velocity, V s , the preshock density of the gas, no, and the 
preshock magnetic field, Bo. The radiation flux is determined by 
the temperature of the star, interpreted as a colour temperature, 
T„, and by the ionization parameter, U. The dust-to-gas ratio, 
d/g is also accounted for, as well as the relative abundances of 
the elements to H. 

A full detailed description of the code is to be presented in 
Contini & Viegas (2007 in preparation). 



Table 2. The radio-mm data references. 



v [GHz] 


Flux [mJy] 


Ref. 


0.843 


10.6+1.2 


Jones (1985) 


1.4 


13.5 


Siodmiak & Tylenda (2001) 


1.6 


15 


Taylor & Seaquist (1984) 


2.7 


33±4 


Purtonetal. (1982) 


4.885 


50 


Aaquist & Kwok (1990) 


5.0 


46±10 


Purtonetal. (1982) 


5.0 


50 


Siodmiak & Tylenda (2001) 


6.2 


65±10 


Purton etal. (1982) 


8.1 


82±7 


Purton etal. (1982) 


8.31 


73.5 


Bastian (1992) 


8.44 


65.3+3.3 


Seaquist etal. (1993) 


8.9 


90±15 


Purton etal. (1982) 


14.5 


90+10 


Purtonetal. (1982) 


14.965 


72 


Aaquist & Kwok (1991) 


90.0 


127±53 


Purton etal. (1982) 


231 


55+16 


Ivison etal. (1995) 
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Fig. 2. The profile of the electron temperature, of the electron density and of the temperature of the grains (a gr -0.2 fim) downstream 
of the reverse (left diagram) and of the expanding shock front (right diagram). The stratification of the most significant ions is shown 
in the bottom panels. For both diagrams, shock front is on the left side. 



Table 3. The infrared data references. 



A [p.m] (phot, band) Flux [Jy] fl Ref 



1.20 (J), 1.64 (H), 2.22 (K), 3.8 (F) 


1.55 - 


2.1-4. 


5-7.1 


Allen (1983) 


8.28, 12.13, 14.65 


7.454 


• 11.93- 


14.67 


MSX6C PSC, Egan et al.* (2003) 


1.235 (J*), 1.662 (H*), 2.159 (K*) 


0.0312 


-0.115 


- 0.594 


Ramos-Farios & Phillips'' (2005) 


0.9 (F), 1.25 (J), 2.22 (K) 


0.0167 


- 0.034 


- 0.741 


DENISE rf 


12, 25, 60 


18.13 


-28.16 


-5.75 


IRAS PSC 


34.6 




16.6 




He et al. (2005/ 


2.5-45.5 SWS 








ISO IDA* 



* 2MASS Photometric System; a: out of Allen, expressed in 10 14 erg s citT 2 A -1 ; b: The Midcourse Space Experiment Point Source Catalog Version 2.3 
(October 2003); c: 2MASS NIR measurements of galactic planetary nebulae; d: The DENIS consortium (third release - Sept. 2005); e: IRAS Catalogue of Point 
Sources, Version 2.0 (IPAC 1986); f: The 35 absorption line towards 1612 MHz masers; g: ISO Data Archive (@ http://www.iso.vilspa.esa.es). 



4. Modelling the spectra 

In this section we present the models which best reproduce the 
whole properties of Hl-36. The models dealing with the ob- 
served line spectra appear in Sect. 4. 1 . They are cross-checked 
by the continuum SED (Sect. 4.2) until a fine tuning of line and 
continuum spectra is found: we refer in particular to the SED in 
the radio range (Sect. 4.2.1) where the turnover frequency and 
the steeper spectral index of the bremsstrahlung at lower fre- 
quencies are highly sensitive to the physical conditions of the 
gas. 

4.1. The line ratios 

The spectra derived from Allen's observations are shown in 
Table |U col. 2; the models calculated for the reverse (m rev ) and 
for the expanding shock (m exp ) are presented in cols. 3 and 4, 
respectively; their weighted sum appears in col. 5. 

The input parameters of the models appear in Table [5] They 
are constrained by the observations, namely, the shock velocity 
V s is roughly determined by the FWHM of the line profiles (in 



agreement with Allen 1983); the preshock density, no, is con- 
strained by the slope of the radio continuum and by the ratios 
of characteristic lines (e.g. [SII] 6717/6731); the preshock mag- 
netic field, Bo, is characteristic of symbiotic systems (Crocker et 
al. 2001); the stellar temperature, T„, is taken from the analysis 
of Allen (1983); the ionization parameter, U, is determined phe- 
nomenologically; the geometrical thickness, D, is constrained 
by the radio continuum slope and by the ratios of the lines to 
HjS; eventually, the relative abundances of the elements are de- 
termined by the line ratios. Unfortunately, in the present case, the 
dust-to-gas ratio (d/g) cannot be directly determined because the 
bump of the dust reprocessed radiation is covered by the emis- 
sion from the summed dust shell and Mira black body fluxes. 

The line ratios calculated by m rev and m exp are summed up 
adopting the weights w shown in Table [5] last row. They are the 
same as those adopted in Fig.[3]to suite the continuum SED. 

The observed line ratios are reproduced by the summed 
model within a factor of 2, except a few lines which are badly 
fitted, e.g. [NI] 5200, since they have a particularly low criti- 
cal density for collisional deexcitation and are very faint at the 
relatively high density of m rev and m exp . We suggest that the 
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observed Hell 1640/ H/3 may be blended with the rather high 
OIII 1640/H/3 ratio, whose calculated value is > 3. An improved 
fit of the [OIII] 5007/[OIII] 4363 ratio could result from a re- 
verse shock with lower densities, which are however incompat- 
ible with the radio SED: as a matter of fact the [Oil] 3727 and 
[OIII] 5007 lines are mostly emitted from the expanding shock. 

Model m rev shows N/H and O/H slightly higher than solar by 
factors of 1.7 and 1.3 respectively, while Fe is strongly depleted 
both in m rev by a factor of 10 and in m exp by a factor of 5: this is 
not surprising as iron, a refractory element, is easily trapped into 
dust grains. 

In Fig. [2] we present the profile of the electron density, elec- 
tron temperature and the fractional abundances of the most sig- 
nificant ions throughout the nebula downstream of the reverse 
(left diagram) and expanding shock (right diagram). It is worth 
noticing that model m rev is matter-bound: this is characteristic of 
the interbinary symbiotic nebulae, where the geometrical thick- 
ness is constricted by the colliding wind region. 



Table 4. The fluxes of the UV and optical spectral lines, normal- 
ized to H/3. 



line 


obs f ' 


m rev 






CIV 1549 


18.78 


19.7 


27.3 


22 


Hell 1640 


7.88 6 


1.1 


5.2 


2.6 


NIII] 1749 


5.15 


2.24 


2. 


2.2 


CIII] 1909 


27.2 


9.3 


18. 


13.0 


[NeV]+OIII 3425 


4.24 


0.13' 


o.i y 


0.22 1 " 


[FeVI] 3662 


0.006 


0.0043 


0.0047 


0.0044 


[OII]+ 3727 


0.576 


0.04 


2.6 


0.99 


[Nelll] 3869+ 


2.12 


1.5 


2.2 


1.78 


[SH]+ 4070 


0.18 


0.38 


0.08 


0.27 


Hy 


0.42 


0.44 


0.45 


0.44 


[OIII] 4363 


0.61 


1.29 


0.80 


1.1 


HeI4471 


0.056 


0.036 


0.023 


0.044 


[Felll] 4658 


0.015 


0.01 


0.015 


0.012 


Hell 4686 


0.55 


0.16 


0.76 


0.4 


H/3 4861 


1 


1 


1 


1 


[ArlV] 4711 


0.02 


0.0013 


0.13 


0.049 


[NeIV]+ 4714 


0.11 


0.14 


0.08 


0.12 


[ArlV] 4740 


0.097 


0.011 


0.15 


0.063 


[OIII[ 5007+ 


19.1 


2. 


35. 


14.2 


[FeVI] 5146 


0.017 


0.007 


0.03 


0.016 


[Fell] 5159 


0.03 


0.039 


0.0017 


0.025 


[FeVI] 5176 


0.03 


0.018 


0.033 


0.023 


[ArIII]5191 


0.005 


0.017 


0.012 


0.015 


[NI] 5199 


0.006 


0.0002 


0.0005 


0.0003 


[FeVI] 5335 


0.011 


0.0018 


0.017 


0.008 


[Nil] 5755 


0.15 


0.43 


0.04 


0.28 


Hel 5876 


0.1 1 


0.18 


0.07 


0.14 


[FeVII] 6087 


0.12 


0.002 


0.005 


0.003 


[OI] 6300+ 


0.4 


0.47 


0.013 


0.3 


[Sill] 6312 


0.08 


0.14 


0.13 


0.136 


[Nil] 6548+ 


1.0 


0.2 


1.4 


0.63 


Ha 


4.5 


3.5 


3.1 


3.4 


[SII] 6717 


0.02 


0.002 


0.057 


0.022 


[SII] 6731 


0.039 


0.0035 


0.10 


0.039 


[Arlll] 7136 


0.24 


0.14 


0.24 


0.14 


[ArlV] 7170 


0.007 


0.023 


0.009 


0.018 


[Oil] 7319+ 


0.55 


2.6 


0.48 


1.8 


[Fell] 8617 


0.0027 


0.0027 


0.0011 


0.002 


[SIII] 9069 


0.33 


0.06 


1.3 


0.5 


Hp abs. rf 




1.13 


0.0075 





a: observed values from Allen (1983); 

b: blended with OIII 1640; 

c: the calculated flux refers only to [NeV]; 

d: H/3 absolute flux calculated at the nebulas, [erg cm -2 s~ ]. 



4.2. The continuum SED 

In Fig. [3] we present the modelling of the continuum. The SED 
is the result of the emitting contributions from the cool and hot 
stars, as well as of the fluxes from the dust shells and of the 
bremsstrahlung from the shocked nebulae downstream of the 
shock fronts, which emit the UV and optical line spectra. 
The bremsstrahlung from the nebulae shows two main peaks. At 
higher frequencies the continuum is emitted from gas collision- 
ally heated by the shock at relatively high temperatures: the peak 
frequency depends therefore on the shock velocity. The peak at ~ 
10 I4 Hz, on the other hand, depends on the volume of gas at tem- 
peratures of ~ 1-3 10 4 K, which is heated and ionized mainly 
by the photoionizing flux: the peak frequency is therefore more 
sensitive to the radiation parameters, T* and U. Also a contribu- 
tion to the SED by synchrotron radiation, not so exotic in objects 
where shocks are at work, can not be ruled out (see Sect. 4.2.1 
and Fig |3]l. 

As a matter of fact, without more constraining data it remains 
puzzling to determine whether the measured flux densities at ~ 
10 n Hz are dominated by synchrotron, by bremsstrahlung or by 
thermal emission from cold dust, also taking into account that 
emission from either mechanism could fluctuate in response to 
variable mass-loss episodes. 

4.2.1. The radio continuum 

In an important paper on radio emission of mass-losing stars, 
Wright & Barlow's (1975) presented an analytic study of the 
spectral flux distribution produced by completely ionized, uni- 
form, spherically symmetric mass loss flow in early type star. 
Though under idealised hypothesis, they tried to interpret the 
radio data for some recently observed objects in the light of 
thermal free-free emission which would result in a power-law 
spectrum S v °c v ' 6 . They also discussed deviations from the 
ideal treatment caused by e.g., non-uniform mass-loss rates or 
by the actual ionization structure of the emitting circumstellar 
envelope. 

In effect, the Purton's radio survey (Purton 1982) pointed out 
that there were some objects (with a high incidence of SS) whose 
radio spectral index was systematically steeper than o'=0.6 (in 



Table 5. The model input parameters. 



parameter 


m rev 




V s (kms- 1 ) 


140 


70 


n (crrT 3 ) 


2.5e5 


3.5e3 


B (gauss) 


le-3 


le-3 


T, (K) 


1.5e5 


1.5e5 


U 


2e-3 


2.5e-3 


D (cm) 


2.8el2 


lel5 


dig 


le-14 


le-14 


He/H 


0.1 


0.1 


C/H 


3.3e-4 


3.3e-4 


N/H 


1.5e-4 


9.1e-5 


O/H 


8.6e-4 


6.6e-4 


Ne/H 


8.3e-5 


8.3e-5 


Mg/H 


2.6e-5 


2.6e-5 


Si/H 


3.3e-5 


3.3e-5 


S/H 


1.6e-5 


1.6e-5 


Ar/H 


6.3e-6 


6.3e-6 


Fe/H 


3.2e-6 


6.2e-6 


log w 


-11.05 


-9.1 
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the range 0.8-1.5), therefore implying an envelope density pro- 
file quite far from the inverse-square one assumed by Wright & 
Barlow (1975). In any event, no attempt was made by Purton to 
fit a detailed physical model to the data. 

Few years later, Seaquist et al. (1984) and Taylor & Seaquist 

(1984) presented a model (hereafter STB model) which ex- 
plained the Hl-36 radio properties in light of the mounting evi- 
dence for the binary nature of symbiotic stars: by assuming that 
the emitting region arises from the portion of the cool compo- 
nent wind photoionized by the WD, they performed a fit to the 
Purton's data with a least-square criterion as a function of the 
so called X parameter (related to the shape of the ionized neb- 
ula), of the turnover frequency v,, and of the flux density at v t . 
The model was successful in reproducing many of the observed 
properties of radio symbiotics, such as the spectral index, the 
spectral turnover, the correlation between the radio emission and 
the spectral type of the cool component. Nonetheless, it was soon 
recognised that the actual radio sources were more complex than 
those assumed in the STB model, and that any derived physical 
property had to be treated with particular caution. Furthermore, 
more recent observations by Seaquist & Taylor (1993) at mil- 
limeter and submillimeter wavelengths showed that several ob- 
jects might differ considerably from the STB model. 

As with regards to the specific case of Hl-36, the STB model 
failed e.g. to reproduce the 843MHz flux presented by Jones 

(1985) . Moreover, the assumed density profile appeared too sim- 
plistic: even the two-layer model discussed by Costa & de Freitas 
Pacheco (1994) referring to the inner structure of the emitting 
nebulae cannot be considered adequate since SS have such im- 
portant density gradients that do not even allow to adopt unique 
values for electron temperature and densities. 
Unfortunately, so far there are no papers which try to interpret 
the radio spectra taking into account the possibility of any out- 
flow from the compact star. 

In the following we fit the existing radio-mm data we have 
presented in Table [3] and Fig. [1] - bottom panel, within the 
colliding- wind scenario depicted in Sect. 3.1. Our approach is 
shown to be self-consistent with the continuum profile of the 
system along the whole electromagnetic range, as well as cross- 
checked by the spectral line emission. 

Three main slopes are seen between 10 s and 10 11 Hz. 
The flattest one, with a spectral index of 0.75, can be explained 
by synchrotron radiation created by the Fermi mechanism at 
the shock front (Bell 1978a,b): a cut-off at the lower limit is 
given by Ginzburg & Syrovatskii (1965) (v ~ 20 N e /B, with 
N e ~10 6 cirT 3 and B=10~ 3 gauss at the shock front), explaining 
the SED profile at ~ 10 n Hz. The unabsorbed bremsstrahlung 
explains the slope of the SED between about 10 ll) and 10 11 Hz. 
At v < 10 10 Hz the bremsstrahlung shows the high slopes which 
are produced by free-free self-absorption. The optical depth t = 
8.24 10 _2 r _L35 v _21 £' (where T is to be measured in K, v in GHz 
and E, the so called emission measure E — J n e n + ds, in cm~ 6 pc - 
Osterbrock 1988) is low for gas at relatively high temperatures: 
however, in the downstream nebulae of Hl-36 we found large 
regions (i.e. large ds) of gas at T ~ 10 4 K (Fig. [2j which make 
r > 1 . The other parameters which vary from slab to slab are the 
density, which is calculated through the compression equation, 
and the geometric thickness of the slab. At this point we would 
like to emphasise that the calculation of the density profile is 
not assumed a priori in our models accounting for shocks, but 
consistently calculated across the whole emitting nebula (Sect. 
3.2). In this case, the thinner the slab the better is the approxi- 
mation of the calculated spectra. Therefore the densities must be 
sufficiently high to give a sensible absorption. The turnover cor- 



responds to different frequencies for different slabs and the final 
bremsstrahlung results from the integrated flux throughout the 
whole length of the nebula. The actual turnover frequency v t in- 
dicates then the transition from the optically thin to the optically 
thick case. 

Eventually, it is worth highlighting that higher densities gen- 
erally refer to the nebula downstream of the shock front fac- 
ing the WD between the stars (model m rev ), while the densities 
downstream of the expanding shock (model m exp ) correspond to 
an optically thin nebula. Therefore, the radio spectrum observed 
in SS may be emitted by the gas between the stars, with an im- 
portant, more increasing contribution from the expanding nebula 
at lower frequencies. Interestingly, this agrees with the variable 
size of the radio nebula, which becomes gradually larger at lower 
frequencies (i.e. 0.6 arcsec at A=2 cm and 5 arcsec at /l=6 cm - 
Taylor 1988). 

log A [/i,m] 

4 2 0-2-4 
-6 — i — i — i — i 1 ' — i — i — i — i — i — i — i — i — i 




10 12 14 16 18 

log v [Hz] 

Fig. 3. The Hl-36 continuum SED, from radio to UV. Thin short- 
dashed line: 250 K dust shell; thick short-dashed line: 800 K dust 
shell; thin long-dashed line: Mira (2500 K) stellar component; 
thick long-dashed line: hot stellar component; dot-dashed line: 
synchrotron; thin solid line: m exp ; thick solid line: m rev (see Table 
|5]for details about the model input parameters); dotted line: the 
jet (see Sect. 4.3). Notice that for m rev , m exp and the jet models 
there are two curves: one is the bremsstrahlung emission from 
radio to X-rays, the other represents the consistently calculated 
reprocessed radiation by dust in the IR. 



4.2.2. The IR continuum 

One of the most intriguing aspect of the symbiotic phenomenon 
in D types pertains to the dusty environment. As a matter of fact, 
these systems show a broad IR excess which has been attributed 
to emission from circumstellar dust since the first IR surveys. 
Early observations already showed that in D type objects the 
dust excesses have colour temperatures near 1000 K (Feast et 
al. 1983); nevertheless today, thanks to the IR astronomy devel- 
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opments, it has been realized that it is difficult to explain the ob- 
served IR spectrum by a single temperature. Theoretical models 
too, confirmed that different "dust" temperatures should be com- 
bined in order to reproduce the NIR-MIR data (e.g. Anandarao 
et al. 1988, Schild et al. 2001, Angeloni et al. 2007c, in prepara- 
tion). 

In Hl-36, the temperature of the dust shells is derived by 
modelling the IR data we collected which, quite surprisingly, 
agree in the overlapping frequency range: this proves that vari- 
ability, even if present as in all symbiotic IR spectra, has not been 
as large as to substantially modify the dust continuum emission 
and then invalidate our results (see Sect. 2, Fig.Q]and Table[3]l. 
It is worth noticing that the cooling rate is very strong at the den- 
sities typical of the region closer to the shock front downstream, 
which corresponds to the maximum temperature of dust: this im- 
plies that the flux from the shell corresponds mainly to the max- 
imum temperature of the grains and can be modelled by a black 
body (bb). As a matter of fact, Fig. [3] shows that the whole IR 
continuum can be well explained by the combination of two bb 
curves, corresponding to temperatures of 800 K and 250 K. By 
comparing the models with the data we obtain the rj factors: they 
depend on the radius of the dust shell, r, and on the distance of 
the system to Earth d (77 = r 2 Id 2 ), being the fluxes calculated at 
the nebula and the data, obviously, measured at Earth. Adopting 
d=4.5 kpc (Allen 1983), we find r= 4.4 10 14 cm and 4.7 10 15 cm 
for the shells at 800 K and 250 K, respectively. Both the derived 
inner shell temperature and radius are in good agreement with 
what found by Allen (1983), who suggested that the cool com- 
ponent is heavily embedded in its own dust shell, specifically the 
hot shell is circumstellar. 

An interesting issue is related to the outer shells: if we relied on 
the Allen's proposed binary separation (3 10 16 cm), then even 
this shell would be circumstellar, namely, well comprised within 
the system. Allen himself admitted such a large binary separa- 
tion was striking, though; moreover, he stated that, translated in 
an angular separation of ~0.5 arcsec according to the suggested 
distance of 4.5 kpc, this binary separation would have been large 
enough to be resolved by careful infrared astrometry, as well as 
by high-resolution optical and radio observations. As a matter of 
fact, few year later Taylor (1988) succeed in resolving the radio 
nebula. Conversely, Bhatt & Sagar (1991) failed in detecting any 
nebulosity in the optical emission lines, arguing that if this was 
actually present, it should be smaller than ~2 arcsec. Corradi et 
al. (1999) finally resolved the optical nebula and found out that 
it appeared also in a continuum image, therefore suggesting that 
it might be a reflection nebula. 

At this point, we would like to stress that the unusual large 
size of the system has been constrained so far by observations of 
gaseous emitting nebulae, and that the binary stellar components 
have never been resolved: this might be a clue that the actual 
binary separation is overestimated, therefore implying that the 
outer shell we found directly by modelling the IR data may even 
be circumbinary, surrounding both the stars. This would not be 
so out of ordinary for D-type SS, since this dichotomy in the 
dust shell distribution (with a hot shell being circumstellar, and 
a cool one circumbinary) has been confirmed by recent studies 
(Angeloni et al. 2007c, in preparation). 

Hl-36 appears as an extreme SS: only further observations 
at higher spatial resolution and in different spectral ranges will 
be able to investigate the actual size and physical properties of 
the dust and gas nebulas. 



4.2.3. The UV continuum 

For sake of completeness we would like to note that longward 
of log v=14.5 the UV data are well fitted by a T,=150.000 K bb 
curve (Fig. [3]), in agreement with the Allen's (1983) derived hot 
star temperature: the WD flux dominates in that range over the 
bremsstrahlung emission. However, it is worth reminding that 
between log v ~ 16.4 and log v ~ 17.7 absorption of X-rays is 
quite strong, therefore implying that X-rays, if present, would be 
soft and mainly due to bremsstrahlung emission arising from the 
high-velocity component we deduce from the IR line features 
(see Sect. 4.3). 

4.3. AjetinH1-36? 

In a previous analysis of ISO SS infrared spectra (Angeloni et 
al. 2007a - hereafter Paper I), we found broad emission lines, in- 
dicating high-velocity components in many objects of the sam- 
ple (Paper I - Fig. 1, Table 4). In Hl-36, the FWHM of these 
line profiles correspond to velocities of 500 kms ; we remind 
that high-velocity components were not observed in the optical- 
UV spectra presented by Allen (1983) which showed FWHM < 
200 kms 1 . Consequently, since the ISO-SWS spectrum shows 
several highly-ionized and broad Ne and O lines, not compati- 
ble with the physical models derived by the optical-UV spectra, 
we decide to include these new results in order to draw out an 
as most large and consistent as possible interpretation. We then 
calculated a further model, characterised by higher shock veloc- 
ities and lower densities, able to reproduce the observed IR lines 
without making the previous fit worse. 

In Paper I we found that, in Hl-36, the [NeVI] 7.65 fim line 
dominates, followed by the [OIV] 25.89 /mi, [Nelll] 15.55 jum, 
[NeV] 24.31 /jm, and [NeV] 14.32 jum lines. 
The model calculated with V s =500 kms -1 , no=10 4 cm" 3 , 
U-S, and D=10 16 cm well reproduces the observed in- 
tensity ratio: in fact we obtain [Nell] 12. 8/[NeVI]7. 65=0, 
[NeV]14.32/[NeVI]7.65=0.47, [NeIII]15.55/[NeVI]7.65=1.3e- 
4, and [NeV]24.31/[NeVI]7.65=0.05, in agreement with the ob- 
served upper limits, respectively [NeV]14.32/[NeVI]7.65<0.5, 
[NeIII]15.55/[NeVI]<0.2 and [NeV]24.31/[NeVI]<0.15. 
Moreover, the same model represents a non negligible con- 
tribution only on the high ionization level lines, particularly 
to CIV 1549, [NeV] 3425, [FeII]6087, and a relatively strong 
contribution to Hell 1640. 

We can now calculate the distance of this high-velocity struc- 
ture from the WD. A WD T„=150.000 K corresponds to a ion- 
izing photon flux of 4 10 26 photons cm" 2 s _1 : this flux is re- 
lated to the ionization parameter and to the gas number den- 
sity, in the radiation dominated zone, by F v (RwdI rj) 2 = Unc. 
Adopting an average density of n=10 5 cm" 3 after compression, 
and a WD radius of 10 9 cm, we obtain log 77 = -16.2. This value 
allows to place the bremsstrahlung emission from this compo- 
nent in the continuum SED diagram of Fig. [3] and to verify that 
its contribution to the radio spectrum is negligible when com- 
pared to the bremsstrahlung from the nebulae. Conversely, on the 
high-energy side of the spectrum, it is interesting to notice that 
V s =500 km s" 1 would correspond to a temperature in the imme- 
diate postshock region of ~3.6 10 6 K, hence indicating that a soft 
X-ray emission may be actually present in Hl-36. 

Eventually, after this consistent analysis, we would like to 
stress the possibility that this high-velocity component in Hl- 
36 is a jet-like feature. This suggestion, if confirmed, would be 
fascinating, as only recently it has been recognising that SS are 
a class of jet-producing objects. Furthermore, the two SS X-ray 
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jets that have been discovered to date (R Aqr, Nichols et al. 2007 
and CH Cyg, Karovska et al. 2007) both show a X-ray structure 
more extended than their radio ones (Sokoloski et al. 2006). This 
would explain why this 10 16 cm jet has not been seen in the 
radio image (with comparable angular resolution) and drive us 
to indicate SS as promising X-ray targets. 

5. Concluding remarks 

In this paper we have analysed the spectra of D-type SS Hl- 
36 within a colliding-wind theoretical framework. After having 
constrained the relative physical conditions by modelling more 
than 40 emission lines from radio to UV, we have been able to 
confidently explain the continuum SED by taking into account 
all the emitting contributions arising from both the stars, the dust 
shells and the gaseous nebulae. A comprehensive model of the 
radio spectra allowed to reproduce the different slopes of the ra- 
dio profile and the turnover frequency, as well as the size of the 
observed nebulae at several frequencies in light of different con- 
tribution from the expanding and reverse shocks. The IR contin- 
uum unveiled the presence of two dust shells with characteristic 
radii and temperatures: the inner shell is confirmed to embed 
the Mira star, while the outer one may be circumbinary, i.e. sur- 
rounds the whole binary system. We believe that the presence 
of multiple dust shells is not a unique characteristic of Hl-36, 
but it may somehow represent a distinctive feature of D-type SS. 
Furthermore, the broad profiles of the IR lines directed us to 
investigate whether a high-velocity component (perhaps an X- 
ray jet) may be present. This insight led us to indicate Hl-36 as 
a promising X-ray target and represents a further support to the 
emerging interpretation of SS as a class of jet-producing objects. 
We then encourage new observations and studies which consis- 
tently take into account the complex nature of SS throughout the 
whole electromagnetic spectrum. 
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